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1. Introduction 

The Wess-Zumino-Novikov-Witten (WZNW) model has received considerable attention 
as an interesting non-rational conformal field theory [|], || and as the Euclidean version of 
AdSs . The 3-dimensional target space of the theory may be parametrized by coordinates 
7 and 7 of its 2-dimensional boundary along with some radial coordinate (p. As for any 
Anti-de Sitter (AdS) geometry, the latter is particularly interesting. Physically, it should 
be regarded as a very close relative of the Liouville direction in 2-dimensional string theory. 

It was long observed |l| that the 3-point functions of the WZNW model coincide 
with those of Liouville field theory ||, |j up to some simple kinematical factors that 
are determined by the SL(2,C) symmetry of the background. In particular, all the 
highly nontrivial (non-perturbative) curvature dependence of the 3-point couplings is 
inherited from the stringy corrections of Liouville theory. A very remarkable generalization 
of this fact was discovered in |7| : Ribault and Teschner determined genus zero correlators 
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for any number N of WZNW primaries in terms of certain 2N — 2 point functions in Li- 
ouville theory. Their proof relies on a relation || between differential equations of both 
theories, namely the Knizhnik-Zamolodchikov (KZ) equations for WZNW models and the 
Belavin-Polyakov-Zamolodchikov (BPZ) equations for Liouville correlators, along with de- 
scent relations provided by the explicit knowledge of the 3-point functions. The relation 
between Liouville theory and the H$ model was further explored in ||, 10, Q |. 



The aim of this paper is twofold. To begin with, we shall re-derive the Ribault- 
Teschner correspondence for genus zero correlation functions using a rather elegant path 
integral computation instead of heavy algebraic manipulations. Our derivation provides 
a completely new view on the map between WZNW and Liouville primaries and on the 
necessity to introduce N — 2 further degenerate fields in the Liouville correlation function. 
The simplicity of our derivation opens the way to various generalizations. Among them 
is the extension of the correspondence to correlators on higher genus surfaces. In fact, 
using essentially the same ideas as for the tree-level computation, we shall derive a precise 
expression for iV-point correlators of WZNW primaries on any closed Riemann surface in 
terms of Liouville field theory. On the Liouville side, the construction involves correlation 
functions with N + 2g — 2 additional insertions of degenerate fields. Further extensions, 
e.g. to the case with world-sheets with boundaries or target space groups of rank r > 1 will 
be briefly discussed at the end of this work. 

Let us now describe the contents of each section and our main results in a bit more 
detail. As we have mentioned, we shall start our analysis by giving a new derivation of the 
Ribault-Teschner correspondence on the sphere. For the convenience of the reader, we also 
include a few comments on the relation between differential equations. Section 3 contains 
the generalization to the torus. For genus g = 1, all the special functions involved in the 
formulation of the correspondence can be easily expressed in terms of Jacobi's 8 function. 
This makes our formulas particularly easy to deal with. In particular, we shall be able to 
demonstrate explicitly how our correspondence intertwines between the relevant differential 
equations. Correlation functions of Liouville theory on the torus satisfy an extension of 



the BPZ differential equation [12| that involves a derivative with respect to the modular 
parameter r of the torus in addition to derivatives with respect to the insertion points of 
fields. On the WZNW side, the story is a bit more complicated. In fact, the usual WZNW 
correlators do not obey KZ-type equations. The problem arises from the zero modes of 
currents which cannot be written as differential operators acting on correlation functions. 



This issue was resolved by Bernard in [13], who had the idea to introduce an additional 
dependence on the choice of some group element. The latter parametrizes possible twisted 
boundary conditions for currents along the /3-cycle of the torus. Not all these dim G 
parameters are actually needed for the formulation of KZ-type equations on the torus. 
The minimal number of required extra parameters is given by the rank of G rather than 
its dimension. In our case this means that one extra twist parameter A is sufficient. Once 
such an extension of WZNW correlators is taken into account, they satisfy the Knizhnik- 
Zamolodchikov-Bernard (KZB) differential equations. By our correspondence, these are 
correctly related to the genus g = 1 BPZ equations. 
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Section 4 contains the main new formula of this work, namely a construction of TV- 
point WZNW correlators on a closed Riemann surface of genus g from correlation functions 
in Liouville field theory, see eq. ( 4.14j ). As we stated before, the latter involves N + 2g — 2 
insertions of Liouville degenerate fields in addition to the N primary fields. It is instructive 
to compare the parameter spaces on both sides of this correspondence. Obviously, the cor- 
relation functions of the WZNW model must be extended by introducing extra parameters 
generalizing the role of the modulus A we discussed at length in the previous paragraph. 
It turns out that on a surface of genus g > 1 we need 2g — 1 new complex coordinates. 
For fixed surface moduli r and insertion points z u , v = 1, . . . , N, an iV-point function of 
WZNW primaries therefore depends on real 3iV target space momenta (j v , p u , p u )u=i,...,N 
and 2g — 1 complex moduli. These add up to 3N + 4g — 2 real moduli. On the Liouville 
side, we count N real target space momenta a u in addition to the (complex) position of 
N + 2g — 2 Liouville degenerate fields and an overall complex pre-factor u. The total 
number of real parameters is therefore 3iV + 4g — 2, just as in the WZNW theory. 



2. Derivation of the correspondence - genus 

Our first task is to explain how the correspondence emerges from a path integral 'definition' 
of the WZNW model. Indeed, we will be able to recover the formula of Ribault and Teschner 
through some formal path integral manipulations. In the second subsection we briefly 
review how the differential equations on both sides of the correspondence are mapped onto 
each other. At genus this is not new, but it will help to understand the corresponding 
analysis at higher genus later on. Finally, we shall comment on the hidden problems of the 
formal path integral approach and the precise interpretation of our results. 

2.1 Path integral derivation of the correspondence 

Let us begin with a little bit of background on the WZNW model. In the most common 
presentation, the action involves the three fields 7, 7 and 4> corresponding to the coordinates 
that parametrize the 2D boundary of and the radial direction, respectively. We will not 
work with this version but pass to a first order formulation which includes two additional 
fields (3 and (3 of conformal weight (hp,0) = (1,0) and (0,hs) = (0,1). Throughout this 
work we shall work in conformal gauge, where the world-sheet metric and curvature are 
determined by some function p through 

ds 2 = \p(z)\ 2 dzdz , y/gK = -4ddln|p| 2 . 

In this particular gauge, the action of the WZNW model takes the following form (see 
e.g. I@) 

S[<j>, 7 , /?] = i- y d 2 w (d<f>d(f> - Hdl -m+^-^K4>- b 2 (3Pe 2b ^ . (2.1) 

Here, the parameter b is related to the level of the WZNW model through b~ 2 = k — 2 and 
the background charge is given by = b. Let us note that the usual action of the 
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model emerges after integration over (3 and (3. The total central charge of the model 
is computed from the level k using 

c(H+) = ^ = 3 + 6b 2 = 2 + (l + 66 2 ) . 

In the last step we split the central charge into the contribution c{(3"f) = 2 from the /?7 
system and the remainder c{4>) = 1 + 66 2 which originates from the bosonic field (f) with 
background charge = b. 

The second ingredients we shall need before we can study correlators of the model 
are the vertex operators of the theory. In the so-called //-basis these read 

Vj{v\z) = \p(z)\ 2A ? M 2J+2 e^W-^eaty+D^) . ( 2 . 2 ) 

In conformal gauge, the factor involving p{z) must be included in order for V to transform 
as a primary of weig 
exponent is given by 



as a primary of weight zero under conformal transformations. The quantity in the 



Af = -6 2 j(j + l) • 

One may consider fi and ft as Euclidean light cone momenta. Similarly, the parameter j 
is related to the momentum in radial 4> direction. In the \x basis, the operator products of 
currents with primary fields are given by 

J a (w) Vj(ji\z) = (w - z)' 1 D a Vj{n\z) for a = ±,0 (2.3) 

where the generators D a of global target space symmetries take the form 

D-=p, D° = -^, D+ = ^ 2 -^±i2 (2.4) 

and similarly for the remaining three generators D a . The fi basis is the most convenient 
one for what we are about to discuss. 

Our aim is to compute the A^-point function of primary fields in the model, i.e. 
> n \ H r N 



Our first step is to integrate out the fields 7 and 7. This is rather easy because they 
appear only linearly in both the action and the exponents of the vertex operators. Hence, 
the integration leads to a simple 5 function constraint on the coefficients of 7 and 7, i.e. 



Af N 



dp(w) = 2ttJ2^ 2 (w - z v ) , d(3{w) = -2vr ^ Jl u 5 2 {w - z u ) . (2.5) 



The distribution 5 2 on the right hand side is normalized such that f d 2 z5 2 (z)f(z,z) = /(0). 
Let us stress that a meromorphic differential j3 with the property ( |2.5D exists if and only 
if the sum ^ jjL v vanishes. Once this condition is met, the integration of 8(3 and 8(3 with 
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respect to the world-sheet coordinate w is immediately performed using the simple auxiliary 
formulas 

0(1/*) = d{l/z) = 2tt5 2 {z) . 

The result is 

N 

/?H = • ( 2 - 6 ) 

A similar equation holds for —(3. The crucial idea now is to re-parametrize (3 using simple 
facts about meromorphic one-differentials on the sphere. To begin with, we recall that for 
any one-differential, the number of poles exceeds the number of zeroes by two. Hence (3 
must have N — 2 zeroes whose locations on the sphere we denote by w = y%. Furthermore, a 
differential is uniquely characterized by the position of its zeroes and poles up to an overall 
factor u. Consequently, we can rewrite (3 in the form 



\U=i {w-Vi ) 

T~riV / 



(3(w) = u " i = 1 ) w yi [ =: uB ( yi ,z u ,w) . (2.7) 



Thereby, we have now replaced the N parameters \i v subject to constraint Ylv flv = ® 
through N — 2 coordinates yi and a global factor u. We can recover the residues [i v of (3 
from the new variables yi and u through 

^ = U U N (z -z ) • ( 8) 

The new variables may be used to rewrite the S function resulting from the integration 
over 7 and j, 

N N 

<5 2 (a/3H-2vr^ / u^ 2 (u;-^)) = 5 2 £>-) 6 2 ((3 - uB (y u z^w)) . (2.9) 
v=l v=\ 

The Jacobian for the transformation from 8(3, d(3 to (3, (3 is trivial on the sphere. Once eq. 
Q2.9| ) has been inserted into our path integral, we perform the integral over the fields (3 and 
(3 to obtain 

( ft V iM*») ) H = l-l^ 2 (E^) / V*e-*'* w (^WW^) x 

N 

x [J \p( z v)\ 2A " |«|- 2 (>+ 1 )|/x l/ | 2 (>+ 1 )e 26 (>+ 1 )^) . 

v=l 

In writing this formula we have also shifted the zero mode of the bosonic field (j) by (j) \— > 
(b — (1/6) In |it|. This removes the u dependence from the interaction term but introduces 
an additional factor \u\ 2 through the coupling of (b to the world-sheet curvature. 

In order to prepare for the second step of our calculation we observe that the expo- 
nential field exp(26c/>) always comes multiplied with |£>o| 2 - For the corresponding term in 
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the action, this has been spelled out explicitly. In the case of the vertex operators, the 
prefactor is hidden in \p JV \ 2 once we insert eq. fl2.8|) . This suggests to introduce a new 
bosonic field <p through 

cp(w,w) := 4>(w, w) + ^ ^2 In \w - yi\ 2 - ^ In \w - z v \ 2 - In |/)(w)| 2 ^ ( 2 -10) 
where the term in brackets is In \Bq\ 2 . Acting with dd gives 



ddip(w,w) = dd4>(w,w) + 



2 



Before we insert our formulas for ip we have to comment on one small subtlety. In fact, at 
various places we will have to evaluate the quantity ln \ w — z v \ at the point w = z u . This 
leads to well-known singularities which need to be regularized. Following Polyakov |l6|], we 
shall use the prescription 

( lim In I it; - z\ 2 ) := -ln|/)(z)| 2 . (2.11) 

\w-*z /fi n 

The right hand side is finite and has the same behavior under conformal transformations 
as the divergent left hand side. Using our change of variables from cf> to (p along with the 
rule ( 2.11 ) it is easy to see that 

|^|2(>+l) e 2fe(>+l)0(^) = y2(3„+l) \p( ZjJ )\- 2 Uu+i) e 26(>+iM«„) _ (2.12) 

After these comments we are prepared to perform our substitution from to (p. We do 
not want to spell out all the details, but let us observe that 



d 2 wd(bdcb = — — / d 2 w 



2?r / T T 2vr 



1 \ 1 ( N ~ 2 N \ 

d(pd(p + —y/gTlip \ - t(Y1 ^(fO -^vM +••• 

' ' V 1=1 v = \ J 



where the . . . stand for terms that do not contain the field (p. There are a few things we 
can read off from this result. To begin with, the background charge of the bosonic 
field <j) gets shifted by an amount AQ = 1/6 to the new background charge Q v = 6 + 1/6. 
Furthermore, the exponents of the vertex operator insertions at z u are all shifted by <p(z u )/b. 
The precise relation is 

Note that the additional power of p is required to match the behavior of both sides un- 
der conformal transformations. More strikingly, new vertex operators exp(— (p(w)/b) are 
inserted at the N — 2 points yj. With a little bit of additional care we can also work out 
the additional (^-independent factors. The result is 



N \ H N 



u=l I v=\ 



II V jM z v ) ) = \e N (u,y u z u )\ 2 6 2 (^^ [ ^ e ^/^(^%^IH 2 ^ ^ 

N N-2 

J] \p(z v )\ 2A " e 2 ( b ^ +1 ^M^) ft \p( yi )\ 2A -Vz» e -i 



x 

v=\ i=l 
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where, following we collected various factors in the function Oat defined by 

N l N-2 N N-2 

e N (u,yi,z u ) = u Yl(z^-Zu)^ 7 YliVi-Vj)^ 7 Yl n^n-ViV^* » ( 2 - 13 ) 

fi<v 8<j ^=1 i=l 

and we introduced the quantities 

= a(Q ¥ , - a) for a = — — or a = a u = b(j ll + l) + — . (2.14) 



Except for the first two factors, the result may now be re-expressed in terms of correlation 
functions of Liouville theory, 

IN \ H N IN N-2 \ L 

( iiVjMzu) \ = 5 2 Q2^)\@N(u, yi ,z u )\ 2 (n^(^) n v -iMvi)) ■ ( 2 - 15 ) 

\u=l I u=l \u=l i=l / 

The Liouville correlator on the right hand side is evaluated with a bulk cosmological con- 
stant hb = 46 2 and the primaries are defined by 

V a (z) = \p(z)\ 2A "e 2 ^ . (2.16) 

Up to an overall constant, our result coincides with the formula found by Ribault and 
Teschner in [^]. 

2.2 From KZ to BPZ differential equations 

The correlation functions on both sides of the correspondence are known to satisfy certain 
differential equations. Throughout this subsection we shall set p = 1. Liouville correlators 
are known to obey the Belavin-Polyakov-Zamolodchikov (BPZ) second order differential 
equations that come with the N — 2 insertions of the degenerate field V-i/2b- m fact, this 
field is well known to possess a singular vector on the second level, or, equivalently, to 
satisfy the differential equation 

d 2 y V_ 1/2b (y) + b- 2 : T(y)V_ 1/2b (y) : = . 

Such an equation holds for each of the N — 2 degenerate fields at the points y{. Using the 
Ward identities for the Virasoro field T(yi), we can convert these singular vector equations 
into N — 2 second order differential equations for the Liouville correlators fl L , 

/ N N-2 \ L 

Dfn L (z u , yi ) = where n L (z u , yi ) = ( f] V au (z u ) ]"[ V_ 1/2b ( yi )) (2.17) 



\U=1 1=1 



and the differential operators were first found by Belavin, Polyakov and Zamolodchikov 
to be of the form 



dy 2 ^i\(yi-^) 2 vi-zv J j£\(vi-vj) 2 y l -yj J 
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Here, di = d/dyi and d v = d/dz v and the conformal dimensions are given by = 
a u (Q — ot u ), as before. 

Let us now turn to the WZNW model. Its correlators are well known to satisfy the 
Knizhnik-Zamolodchikov (KZ) differential equations. These emerge from insertions of the 
Sugawara singular vector 

T(w) +b 2 (^: J°(w)J°(w) : ~( : J-(w)J + (w) : + : J + {w)J-{w) : ) \ = 

into WZNW correlation functions. Since we can insert this at any point on the sphere, 
we certainly get an infinite number of equations. But these are not all independent. In 
fact, evaluation of the residues of the first order poles at w = z v gives ./V independent 
first order equations. Hence, there are two more equations on the WZNW side than on 
the Liouville side. Closer inspection shows that the N — 2 BPZ differential equations 
correspond inserting the Sugawara tensor at the N — 2 points w = y% on the world-sheet of 
the WZNW model. At these special points, the differential equations from the Sugawara 
singular vector simplify considerably. To see this we need the operator product expansions 
between the currents and J° and the vertex operators in the \i basis (p.3|). From the 
current Ward identities and the very definition of the points yi we conclude that 

N 

By similar reasoning, the current J°(yi) is easily seen to simplify as follows 

jOfyA = _y^_^ ^_ =_A 

^[Vi- z v dfi u dyi 

where the second equality follows from the explicit chance of variables relating fj, v with y, 
and u. Once we insert all these results into the Sugawara singular vector, we obtain the 
following N — 2 differential equations 

Dfn H (z v ,^) = where Q H (z u ,fi u ) = / JJ V jv {fJL v \z v ) ^ . (2.19) 



Using the relation (2.8), the differential operators D^ 1 for H% can be shown to acquire the 



following form 

where AjJ = —b 2 j u {j u + 1) are the conformal dimensions of our WZNW primary fields, 
as before. Note that does not depend on any of the yj with j ^ i. In this sense, 



the transformation ( |2.8| ) leads to a separation of variables [17]. Let us stress that the 



derivatives d u in eq. ( |2,19j ) are to be taken with fixed \i v whereas the variables yi are kept 
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fixed when evaluating d v in the BPZ differential equations above. Whenever it is relevant, 
we shall explicitly distinguish between the two derivatives d u , 



dz v J " \dz u 



fj, y 
It is certainly possible to express e.g. all the derivatives in terms of <9^,<9j and d/du. 
We shall only need the special combination 

Si ■= = E^tW + *)-E — (*-^) • ( 2 - 21 ) 

With this auxiliary formula it is then straightforward to check that 

(e^Dfe N - Df) n L (z„, yi ) = o. (2.22) 

The correspondence between differential equation was an important ingredient in M for 



proving the relation ( 2.15 ). We have derived the latter within the path integral approach 



and therefore the differential equations are guaranteed to be mapped onto each other. 
2.3 Comments on the path integral derivation 

It is well known that the path integral definition of both the WZNW model and of Liouville 
theory has some issues that are related to the non-compactness of the background (see e.g. 
Jj~8| ] for a review). If one splits all fields into their zero modes and fluctuations, one can 
integrate out the zero modes. The procedure results in expression for the correlators which 
are either divergent or hard to give an exact definition. In the case of Liouville theory, for 
instance, one obtains 

(II e 2a ^~ z ^) L = f V<fe~ S ^ n e 2 ^- 2 ')^-^ (p B J d 2 we 2b ^j 

where sb = — a v ^ n ^ ms formula, the integration over the fluctuation field (p 

is weighted with the measure of the free linear dilaton theory, but the integrand contains 
additional insertions of screening charges. The latter are raised to some possibly non- 
integer power s and they are multiplied with coefficients that diverge for positive integer 
powers. These features have a simple explanation. In our non-compact target space, the 
correlation functions are expected to possess poles in momentum space which come from 
the integration over the infinite region of the target space where the interaction is negligible. 
The path integral computation we have just sketched detects (some of) these divergencies 
and can be turned into a rigorous computation of the associated residues. 

If we are not willing to give the path integral any more credit, then our results above 
only imply that 

IN V 07 N 

(HVjMz^S;) = 6 2 (J2^)\ @ N(u, yi ,z v )\ 2 x (2.23) 

\ v=l I LD<t> v=l 

I N N-2 \ 

x ( II V av {z v ) J] V„ 1/2b (y t ) S' v j 
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Here, and 5* denote the screening charges of the WZNW model and Liouville field 
theory, respectively. They are given by the following expressions 

S s ^ = - J d 2 w f3(w)P{w)e 2b(l)iw ' 1J,) , S' v = fd 2 we 2b ^ w '^. (2.24) 

The correlation functions on both sides of the equality are to be computed in the free linear 
dilaton theory. On the left hand side, we use the ^7 system with central charge c = 2 and 
a linear dilaton with background charge = 6. On the right hand side, the correlator is 
to be computed for a linear dilaton background with Q v = 6 + 1/6. 

Correlation functions in the WZNW model and in Liouville field theory are well known 
to possess a second series of poles that are not explained by insertions of the screening 



charges ( 2.24 ). The residues of these poles can still be computed from free field theory 



with the help of so-called dual screening charges. For our two models these read 

d 2 w p{w)^p{w)V , = jd 2 we^ w ^ . (2.25) 

Note that the exponents 1/6 2 = k — 2 can be integer for integer level k > 2. We have indeed 
checked by explicit computation that equation ( [2.23 ) remains valid if the screening charges 



and £7* are replaced by the dual ones. Similar calculations for the H% model involving 
both screening charges and can be found, e.g., in []l9| (see also references therein for 
earlier works on this subject). We shall perform free field computations for correlators on 
surfaces of genus g at the end of section 4. 



3. Generalizing the correspondence to the torus 

Motivated by our rather simple path integral derivation, we would now like to extend the 
-ff^-Liouville relation beyond tree level to higher genus correlators. We shall study the 
general case in the next section and restrict ourselves to the torus for now since many of 
the formulas can be made very explicit at g = 1. 

3.1 The H3 model on the torus 

We start from the WZNW model with level k and compute TV-point functions on a 
torus with moduli parameter r using the first order formulation in terms of <j), (3, 7. As 
explained in the introduction, we would like to introduce a bit more freedom by admitting 
some non-trivial boundary conditions for (ft and the (5^ system. To be precise, we assume 
the fields (f>, 7 and (5 to satisfy 

(3{w + m + nT)= e 2mnX f3(w) , 7 (w + m + nr) = e~ 2mnX ^(w) , (3.1) 

27rnImA 

4>(w + m + nr,w + m + nr) = <p(w) H (3.2) 

6 

with n,m G Z and A some complex parameter. Once such twisted boundary conditions 
have been introduced for the field /3, the conditions on 7 and follow if we require the 
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action to be single valued. Of course, we assume similar twisted boundary conditions with 
twist parameter A to hold for the anti-holomorphic components of the f3j system. 

There is one term in the action, namely the coupling of the field cf> to the world-sheet 
curvature, that requires a bit of additional care. Since our field (j) is multivalued, the term 
^JglZcj) cannot possibly give the right prescription. Instead, we must decompose <j) into 
a twisted zero mode part <f) so \ (the index 'sol' stands for solitonic) and a doubly periodic 
fluctuation cp q , i.e. 

2ir 

4> = 4>q + 4>&o\ where 4> so i(w,w) = — -ImAIrmu , (3-3) 

where, by construction, <p q now satisfies 4> q {w + m + nr, w + m + nf) = (f) q (w). The linear 
dilaton term in the action couples the world-sheet curvature to the single valued fluctuation 
field (j) q rather than to (j) itself, 

SldM = / d 2 w^ll<P q . (3.4) 
All other terms in the action (^]l]) remain the same as before. Similarly, the expression for 



vertex operators (2.2) does not require any modification. 

To a large extend the path integral computation of the iV-point correlation function 
follows the same steps as before. The quantities we would like to compute are given by 

IN \ H r N 

( II V M*») ) = / ^ A ^ A 7^ A /3 e" 5 ^] J] VjMzu) , 

\ u=l I ( AjT ) J v=\ 

where the integration is to be performed over all field configurations on a torus with 
modulus r satisfying the boundary conditions ( |3.1[ ) and ( |3.2|) stated above. Note that 
our correlation function is not yet normalized by dividing the partition function Z H . We 
shall further comment on this below. After integration over 7 and 7 we obtain the same 
condition ( |2.5| ) for the derivatives of (3 and (3 on the sphere. But this time it has different 
consequences for (3 and f3. In fact, we need a bit of preparation before we are able to spell 
out the analogue of the important equation fl2.7|) . 

On the torus, the integration of equation ( p.5|) will lead to a new function that can be 
constructed out of Jacobi's theta function 

9(z) = e i7r ( n+ i )2r+27ri(n+ 5)( 2+ 2) , 9{z + m + nr) = -e^ m ^ z+m ^e{z) . (3.5) 

Out of 6 we can build a new function <j\ with a simple pole and the same periodicity 
properties that we required for j3, 

9(X-(z-w))9'(0) 

a ^ w) = 9(z-w)9(X) ■ (3 - 6) 

Indeed, from the shift properties of Jacobi's theta function 9 it is easy to derive the following 
property of a\, 

a x (z + m + nT,w) = e 27rinX a x (z,w) . (3.7) 
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Now let us return to the integration of the two equations fl2.q ) . The right hand side tells us 
that the twisted meromorphic differential (3(w) possesses N poles in the positions w = z u 
with residues (i u . The solution to these conditions is unique, as long as the twist parameter 
A does not vanish. It can be written down in terms of the <j\ as 



P(w) = y2vucr\(w,z u ) = u 'tr 1 — =: uBi(yi, z v ;w) ■ (3.8) 

u=l Uu=l d ( W 



The second equality is a torus version of Sklyanin's separation of variables and it requires 
a few comments. On a torus, meromorphic one-differentials possess the same number N of 
poles and zeroes. Moreover, the positions w = z u and w = yi, i = 1, . . . , N, of both zeroes 
and poles determine the differential up to an overall factor u. For our correspondence 
between the WZNW model and Liouville theory, it is again crucial to parametrize p(3 
through u and yi rather than A and fi u . The relation between the two sets of parameters 



can be worked out easily (see also [20. 21]) 



What we have shown so far is that the integration over the fields 7 and 7 in the WZNW 
model leads to the following 8 function 

N 

S 2 (d(3(w)-27rY,^S 2 (w-z u )) = \detd\fS 2 (P(w)-uB 1 (y u z u ;w)) . (3.10) 

-2 



This replaces the related formula ( |2.9| ) in the genus zero analysis. The factor | det<9|^ is 
the Jacobian that arises when we change from 5 2 (d(3 • • • ) to 5 2 (f3{w) • • •)■ We have placed 
a subscript A in the Jacobian to remind us that d is considered as an operator on twisted 
one-differentials. Let us now perform the integration over (3 and (3 to obtain 



N ' H 



,,-1 / /, I det9| A J 



X 



u=1 1 (A,r) 

N 



x ni^)i 2A 'i n i~ 2( ^ +i) i^i 2( ^ +i)e2fe( ^ +i) ^ ) • 

As in our genus zero analysis we have shifted the zero mode of the field <j) to remove the 
\u\ 2 from the interaction term. Because the Euler characteristics of the torus vanishes, the 
path integral is multiplied with a factor |it|° = 1. 

We have now reached the point at which we change variables for the remaining in- 
tegration over cj) such that the highly non-trivial factor \Bi\ 2 also gets removed from the 
interaction. In complete analogy to the genus zero case we introduce 

tp(w,w) := ^,^) + i(^ln|^- yi )| 2 -^m|^-^)| 2 -ln|^)| 2 ) • (3.11) 
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By construction, the new field (p is periodic under shift by n + mr, even though the original 
field (j) an d the 6 functions are not. It will be advantageous to replace \8\ 2 functions by 
some function F, 



F(z,w) ■= e T 2 y y " 



6{z-w)\ 2 



no) ■ (3 - 12) 

which is easily seen to be invariant under shifts by integers and integer multiples of the 
modulus r, 

F(z + m + nT, w) = F(z,w) 

The new function F allows us to rewrite the relation between <p and (ft in terms of the single 
valued fluctuation field 4> q = (f) — <j) so \ that we introduced in eq. ( |3.3| ) , 



<p{w,w) :=<f> q (w,w) + ^ ^hiF(w i y i )-^hiF(w,z v )-hi\p(w)\ 2 + Si\ (3.13) 

with Si = 7^X4 2/f ~~ J2u z i)- Let us also observe that the decomposition of the field (j> 
into 4> so \ and <p q is such that their contributions to the kinetic term decouple, 

J d 2 wd(pd(i)= J d 2 wd4> so \d(i) so \ - ^- J d 2 wd4> q d(t) q (3.14) 

with ^- J d 2 wd(j) so id(f) so i = . (3.15) 

Now we can proceed exactly as before, with the function F replacing \z — w\ 2 . In particular, 
the properties of F can be used to evaluate dd(p as 

dd<p(w,w) :=ddcl> q (w,w) + ~ (^5 2 (w- yi )-^5 2 (w-z^-^dd\n\p(w)\ 2 , (3.16) 

a result that agrees exactly with the corresponding formula at genus zero. The outcome of 
a short and straightforward computation is 

/ n VM*>) ) H = e-^ |e rVy )|2 r^-bs^^^v^w^) 

\u=l / (AiT) l det ^A J 

N N 

x II \p(z u )\ 2A »e 2 ( b ^ +1 ^M^ H \ p ( yi )f*-xm e -frM 

u=l i=l 

where the pre- factor G^v is given by 

N N N 

\Q 9 =\y u z v , T )\ 2 = WF{z^z u )^WF(y uyj )^ Yl F( Zfl , yi y^ . (3.17) 

\i<v i<j (J.,i=l 

So far, the correlators on both sides of the equations are not normalized. But using the 
following formulas (see \22\ and our discussion before eq. ( [4.13 ) later on) for the genus one 
partition functions Z H and Z L of the WZNW model and Liouville theory, respectively, 

7r(ImA) 2 

g b 2 T 2 I 

ZH = VWW ' ZL = V^\v(r)\ 2 ' (3 - 18) 
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relation between the level k = b 2 + 2 and the background charge = b + 1/b and 
the same bulk cosmological constant Hb = 46 2 . For each primary field Vj in the WZNW 



along with the formula det(9\ = 9(X)/r](q) for the determinant of d on A-twisted dif- 
ferentials, we can recast the relation between correlation functions in both theories in a 
particularly simple form 

/ N \ H _ IN N \ L 

zh{ n vm*v) ) = le^w^i 2 -^/ n v av { Zv ) n v_ 1/2b ( yi ) \ . 

\ v=\ I ( A) r) \ f=l i=l / 

(3.19) 

Our genus one relation is similar to the tree level result. Once more we managed to 
compute all correlators of WZNW primaries in terms of Liouville theory with the same 

h 

model there appears one Liouville vertex operator V a where a = b(j + 1) + 1/26, as on the 
sphere. Additionally we have to insert degenerate Liouville fields, but now we need two 
more than on the sphere, i.e. there are N degenerate fields inserted. Their positions are 
determined through the light cone momenta fi u of the WZNW vertex operators and the 
twist parameter A. 

The attentive reader might be a bit surprised not to see any factor implementing the 
conservation of [i momentum, as on the sphere. Its absence is directly related to the fact 
that the correspondence has been worked out for non-zero twist parameter A. One would 
certainly expect to recover \i momentum conservation in the limit A — > 0. In order to see 
how this works, let us go back to relation ( |3.8[ ) and insert the expansion 



This leads to the formula 



i e'(w-z) 



= ^^+^ ( vft„lll0|>-2„) + ... . 

V V 

For (3 to stay finite in the limit A — > 0, we need to assume that the total [i- momentum 
Y2 V l*>v tends to zero when we send A — > 0. In fact, if Yl v vanishes fast enough, we obtain 
well defined expressions at A = 0. 

3.2 Relation between differential equations on torus 

Having established a simple relation between correlation functions of the WZNW model and 
Liouville field theory on the torus it seems worthwhile to look once more at the differential 
equations that determine correlators in both models and to check that our relation ( |3.19| ) 
correctly intertwines between them. 

Let us start on the side of Liouville field theory. The vertex operator V_i/2b belongs to 
a degenerate representation with a null vector (6 2 (L_i) 2 + L__2)| — 1 /26) on the second level. 
Since we have such degenerate fields in our Liouville correlation function, we obtain N 
second order differential equations, 

in N \ L 

Df(n L )(z„, yj ,T)= 0, ^ L {z v ,y U T)= {T\V av {z u )\\V_ l/2b { yj )\ . (3.21) 

\u=l j=l I 
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The differential operators D\ for surfaces of higher genus were worked out by Eguchi and 
Ooguri @, 

D > = " 2 $ + A -i 2 '" + g(« ! '- ! '' ) i|" A -i a ' ?<! " ,w) ) + 



where the special function £ and the constant r/i are constructed from Jacobi's function 
through 

9'(z-w) 16"' (0) 



0(z- w) ' 771 6 #'(0) ' 

Let us now address the differential equations obeyed by the correlation functions of the 
WZNW model on the torus, which were first worked out by Bernard and are known as 
Knizhnik-Zamolodchikov-Bernard (KZB) equations. These equations are obtained by in- 
serting the Sugawara singular vector T(w) — b 2 : J a J a '■ (w) = into the WZNW Appoint 
correlation function Cl H (see eq. ( |3.29| ) below for the definition of Q ). Ward identities for 
currents and the Virasoro field then give |0|, 23] 



(3.23) 



N , „ \ f) 

b 2 s(w) K(i/a,^)^- - *fAZ(vi>*»)) +27ri fr 



nn H = o 



(3.24) 



Here, we introduced n = |#(A)| 2 and 



S(w) 



d_ 



J°(w) 



J-(w)J + (w) + J + (w)J-(w)) , 



N 



N 



J T (w) = ^2 a±x(w,z u )D* , J°H =J2t(v,*v)Dl 



(3.25) 



(3.26) 



and the differential operators D^j and are the same as in eq. ( p.4| ). Let us briefly recall 
the reason why the KZB equations contain a derivative with respect to the twist parameter 
A. These terms arise from the Ward identities of currents. In fact, it has already been 
observed by Eguchi and Ooguri in |T^] that the insertion of the zero modes of currents into 
correlators cannot be converted into differential operators acting on the usual untwisted 
correlation functions. It was Bernard's idea to fix this problem by introducing a dependence 
of conformal blocks on additional parameters. On the torus, he suggested to insert a group 
element g into the trace. This has the effect of twisting the boundary conditions for currents 
under shifts by multiples of r. Our boundary conditions correspond to the special choice 
g = exp(— 2iri\jQ ). Actually, it had been observed by Bernard already that a single twist 
parameter A suffices on the torus. 

As in subsection 2.2, our strategy now is to evaluate the KZB equations at the iV 
special points y» and then to compare the result with the iV differential equations for the 
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Liouville correlator. From the relation ( |3.9D between //j,A and yi,u one may derive 

These relations between derivatives can be inserted into our formulas for the differential 
operators J (w) and J°(w), evaluated at the points w = yi, to obtain 

N £ d d d 

When we plug these expressions into S(yi) we find 

Sfa) = i , (3.28) 
dyf 

just as on the sphere. In conclusion we have shown that the KZB equations lead to the 
following N differential equations for the WZNW iV-point functions Q H 

Dfun H = o , n H { Zu ^ u ,\) = (jlVjMzin , (3.29) 

\»=1 / (A,r) 

A H = ^+E(^^)^--A2^(w^))+27ri|; • (3.30) 

Let us recall that the derivatives d/dz u = <9j, in are still taken while keeping A and \ii 
fixed, in spite of the explicit appearance of derivatives with respect yi. 

In order to verify consistency of the two sets of equations with the proposed relation 
between Liouville and WZNW correlation functions, we rewrite the latter in the form 

un H = \r)(r)\ 2 \e' N \ 2 n L , (3.31) 

N N N 

e' N = e'(o)£ Yle^-z^Heiyi-yj)^ J] 9{z v - yi )-^ . (3.32) 

v<[i i<j v,i=l 

In comparison to the earlier version, we have absorbed a factor into a re-definition of |@aH 2 
and then expressed the new |0^| 2 in terms of \6\ 2 rather than F. Our result on the relation 
between correlation functions therefore implies 

(r^r^Mv) = A L ■ (3-33) 

This can be checked indeed by a lengthy but straightforward computation. In the process 
it is important to replace all the derivatives d v = in the differential operators 
through derivatives d u = where the latter are taken while keeping y, and u fixed. More 
precisely, we replace ^ = ZiVh z v)d„ in Df by 



'V 

N 



Si = tin, *>) {du + d i)~Yl t(vi> yi) & - 9 j) . ( 3 - 34 ) 
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where di denote differentiation with respect to yi, as before. It is easy to verify that 5i 
satisfy 5i/j, u (yi, z u ,u) = and 5i\(yi, z u ,u) = 0. Once this issue is taken care of, we can 



show that the differential operators indeed satisfy ( 3.33 ). The details of the computation 
are presented in appendix |A[ 

4. Generalization to arbitrary genus 

Equipped with the experience from genus g = and g = 1, we now address the WZNW 
model on an arbitrary closed surface of genus g. Most of the analysis follows the ideas of 
previous sections, but the details require considerably more background concerning differ- 
entials on higher genus surfaces. We provide the most relevant details in Appendix [B]. Since 
this section contains our main new result, we shall derive it first through our path integral 
arguments and then verify with the help of free field computations where both sides of the 
proposed relation possess the same residues. Comments on the relation between the KZB 
and BPZ differential equations are deferred to the concluding section. 

4.1 Path integral derivation 

Suppose we are given some compact Riemann surface £ with moduli parametrized by the 
period matrix r = (r^-). Since our fields are going to be multivalued as in the torus case, 
it is appropriate to pass to the universal cover £ of the surface right away. Using the 
famous Abel map, we embed £ into C 9 . From now on, we shall think of our fields as 
being defined on the image of the Abel map and hence consider them as functions of g 
complex coordinates w k , k = 1, . . . ,g. In close analogy to the genus g = 1 case, we allow 
for nontrivial twists along the /3-cycles, i.e. 

0(w k + r kl n l + m k \r) = e WA </3(™ fe |r) , 

l{w k + T kl n l + m k \r) = e- WAi 7 (u; fc |r) , (4.1) 



/ i\ if i \ 27rn i ImA/ 

(f){w k +T k in +m k \r) = (p(w k \r) H . 

b 

The complex parameter A^ represents the twist along the /3-cycle (5 k . Thereby, we have 
introduced g complex parameters. 

As in our discussion of the theory on the torus, spelling out the coupling of </> to the 
world-sheet curvature requires to split </> into a twisted zero mode </> so i and a single valued 
fluctuation cp q . The twisted zero mode 4> so \ is now given by 



^sol 



2-77- 2tt f' w 

— Im\ k (Imr)^ 1 bW = — ImA fe (Imr) fe / Im / J (4.2 



where u l , I = 1 ,g, is a basis of holomorphic one- forms and indices I, k are raised and 

lowered with the trivial metric. The linear dilaton term couples the world-sheet curvature 
1Z to the doubly periodic fluctuation field <p q = <f> — (j) so \ in the same way as on the torus, 
see eq. (3.4). 



-17- 



So far, setting up the path integral for WZNW correlators on a surface of genus g was 
a straightforward extension of the torus case. But there is one important modification. As 
is well known, the A-twisted differentials p(3 and p/3 possess g — 1 zero modes. These give 
rise to g — 1 additional moduli if we decide to fix the value of the p(3 and p(3 zero modes 
and to extend our path integral only over the remaining fluctuations. Our aim therefore is 
to compute the following iV-point correlation functions 

in \ H r N 

( II VjM*») ) = / ^V^ A 7^ A /? e~ s ^ J] VMz„) 

\ v=\ I (A,tJ7,r) "=1 

over a Riemann surface with genus g. The symbol reminds us not to integrate over j3 
zero modes. We parametrize the latter by g — 1 coordinates w = (za a ,cr = 1, . . . , g — 1) and 
place an explicit subscript on the correlator. The physical correlation functions may be 
recovered in principle through a finite dimensional integral over w a . 

Integration over 7 leads to exactly the same expression ( |2.5| ) for the derivative of (3 
as on the sphere and torus. But the corresponding (3 takes a different form. From the 
knowledge of its derivative and the boundary conditions (4.1) we may conclude that (3 
must have the form 

N g-l 

P( w ) = ^p>u<?\{w,z u ) + ^2,w a u^{w) (4.3) 

v=l o=l 

where ui^. denote a basis of A twisted holomorphic differentials and the function <J\(w, z) is 
the following differential, 

^ - m ej y^ • 

On the right hand side, we can use any odd spin structure 5. The 6 function 0g and the 1/2- 
differential h$ are defined in appendix ||. Using properties of these objects it is possible to 
show that a\(z,w) has a simple pole at z = w with residue Res z=w a\(z,w) = 1. Moreover 
the differential <j\ satisfies the same periodic boundary condition as f3(w). Before we go 
on, let us point out that we had to fix the (3 zero mode in order to be able to reconstruct 
(3 from its derivative. Explicit formulas for the g — 1 twisted holomorphic differentials can 
be found e.g. in [^] . 

The rest of our analysis proceeds essentially as before. A meromorphic one-differential 
with iV poles is known to possess ^ + 2(^ — 1) zeroes. Hence, the analogue of the separation 
of variables formula (|2.7|) for a surface of genus g > 1 is given by 



«-) = „ n£g^ii^ =! uBAktVtW) (45) 

with yi parametrizing the zeroes of (3. It seems that this formula has not appeared in the 
literature before. In order to reproduce the correct boundary conditions ( f4.ip , yi have to 
satisfy the condition 

A/ = £ / m - / - 2 / w, , (4.6) 

k=i J™ u=l^ w J(g-l)w 



- 18 - 



which corresponds to the condition ( |3,9[) for the genus one case. Notice that this condition 
does not depend on w. The function a(w) in eq. ( |4.5| ) is a g/2-differential without zeros 
nor poles. Its definition is reviewed is eq. ( |B.13| ) of appendix [B] along with the construction 
of the prime form E (see eq. ( |B.7D ) and the Riemann class A. The function a(w) is needed 
in order for the right hand side of (^1]) to become a one-differential with the correct zeros 
and poles. 

There is another point we would like to stress. Note that the equations ( |4.6| ) impose g 
constraints on the position of the zeroes j/j. Thereby, they define a N + g — 2-dimensional 
hyper-surface in the configuration space of y^. These hyper-surfaces sweep out the entire 
configuration space as we vary the g twist parameters. Changing the zero mode parameters 
w a while keeping A fixed also moves the position of zeroes yi, but this motion takes place 



within the the hyper-surface defined by A since eq. (4.6) is independent of w. Hence, the 
infinitesimal changes of the twist parameters A and the zero mode parameters w span a 2g— 
1 dimensional subspace of vectors tangent to the configuration space of j/j. Together with 
the shifts of the N light-cone momenta fi u , we thereby generate independent moves of all 
the zeroes m and of u. Let us also note that the light-cone momenta may be reconstructed 
from yi according to 

= u J . (4.7) 



After these comments we can continue with our computation of WZNW correlators. Once 
the trivial integration over non-zero modes of (5 and (3 has been performed, we re-define 
the bosonic field, 

tp(w,w) := 

/iV+2(9-l) N 
<j>(w,w) + -[ YI ^\E(w,y i )\ 2 -Y l n\E(w,z v )\ 2 + 2ln\a(w)\ 2 -ln\p( W y 2 

y i=i u=i 

As on the torus, we may replace the multi-valued prime form E and a by doubly periodic 
functions and thereby rewrite ip in terms of the single valued fluctuation field (ftq — (f? — ^soi? 

1 /iV+2( fl -l) N 

<j> q (w,w) + — \ Y ^F(w iyi )-Y lnF ( w ^^ + 2lnH ( w )- ln \p( w )\ 2 + S 9 



i=l v=\ 



The functions F(z),H(z) are defined in appendix [b| and the constant S g is a shorthand for 
the following expression 

rvi ryi rzu ryv 

S g = 2irY lm u l (Imr)^ 1 Im / oj k - 2^Y Im / ^ ( Imr )z~fc 1 Im / ujk ~ 

i J WO Jw Q v Jw Jwq 

-47rlm / uj 1 (Imr)^ 1 Im / 

J(g-l)wo J(g—l)wo 
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which is independent of wq. If we act with dd on <p we obtain the same expression ( 3.1 6| ) 
as on the torus. This uses that the contributions from the non-holomorphic part cancel 
each other. Furthermore, ddlna(z) = because a has neither zeros nor poles. 

After inserting the shift of variables from <j> to <p, we can simplify the resulting ex- 
pressions pretty much in the same way as for the torus case. Once more, the kinetic term 
for the field <f> splits into the sum (3.14) of a constant term and a kinetic term for the 
fluctuation field cp q . The former is given by 

Ssoi = ^J d2w ^sol^soi = JlmA' (Imr)-, 1 ImA fc . (4.8) 

A second auxiliary result concerns the contribution from the linear dilaton term. After 
the change of variables it provides us with a linear dilaton term for tp and the following 
additional terms, 

/JV+2( S -1) N 

8 d2w zT\ 2 ^F{w, yi )-Y,^F{w,z u ) + 2\nH{w)-\n\p{w)\ 2 + S g 

J \ 1=1 v=\ 

N+2{g-l) N 
= Yl H\P^)r 1 H(y i ))-^2H\p(z u )\- 1 H(z u )) + (l-g)S g + -U g , 

i=l u=l 

where (see [pif ) 

U g = ^2 / d 2 wd 2 y^gJ^n(w) V r g^TZ(y)lnF(w,y) . (4.9) 

All these expressions can be verified using formulas from Appendix |B|. Collecting all the 
above facts, our correlation function is given by 

IN \ H IN N+2(g-l) \ L 

U,yi, Zy^T 

\ v=\ I (\,vj,t) \ v=l i=l I T 

where the pre- factor Otv takes the form 

N N+2(g-l) 

\%(u, yi , Zu ,T)\ 2 = duivo^n^)" 1- ^ n H (yk) i+ ^ >< (4-io) 

i=l k=l 
N JV+2(ff-l) N N+2(g-l) 

x\\F{z r ,z s )^ JJ F(y k , yi )^\\ J] F{z r ,y k y^ , 

r<s fc<Z r=l fe=l 

and we collected all the remaining terms in the quantity C, 



,2+3^ 9 |det5| 



A 



Before we conclude, let us observe that the constant C may be written as the ratio between 
partition function Zq = Zj^ of a linear dilation with background charge Q v and the 
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product Zq = \Z^\ 2 Zj^ where Z^ is the partition function of a chiral (3j system. In 
order to see this, we recall that the partition function in both model and Liouville field 
theory acquires it's leading (divergent) contribution from the asymptotic region where the 
interaction is negligible. Hence, we have 



= |det<9|- 2 e-^' j V(f)qe -^S<Pu,{84, q d4> q +^V3K4> q ) (4n) 
To further re- write the partition function Zff , we must now change the background charge 



from Qrp to Q v . We can achieve this using a result of [24] on correlation functions in a 
linear dilaton theory. When applied to the path integral in the formula for Zq , it formally 
reads 



/ 



V(j) q e 



J_ r,r- „■(, ■),,;,.,. . /tit? e < f 4 f/ '< 

\detd\ 2 



jv^-*!*^ 6 *** 3 ?**-*) = e^) Ua Z^ (4.12) 



where in the process of the calculation we inserted the explicit expression for the background 
charges and Q v . Combining the previous two equations, we have shown, 

Z* = \detd\fe~ Ssol+ ^ + ^ U3 Z{; . (4.13) 

In conclusion, our final result for the relation between normalized correlation functions in 
the iig model and in Liouville field theory reads 

1 I N \ H I N N+fy-l) \ L 

Vh ( EI VjM*>) ) = \® 9 N (u,y u z u ,T)\ 2 -j. ( J] V ai/ (z u ) J] V-iMVi) ) 

\ u=l I (A>W)T ) Z \ u=1 i=1 I T 

(4.14) 

where the function Q 9 N is given by eq. ( |4.10 ). Let us also recall that the correlation functions 



of primaries Q2.2p in the WZNW model at level k = b 2 + 2 depend on the g twists and 
on g — 1 zero modes w\ in addition to the surface moduli. On the Liouville side, we compute 



the correlation functions of primaries (2.16) with a background charge Q = Q v = b + l/b 



and with bulk cosmological constant hb = 46 2 . The momenta a v are related to j v through 
eq. ( p. 14 ) . The remaining momenta fi u in the WZNW model along with the 2g — 1 moduli 



Afc and w\ determine the insertion points yi of N + 2(g — 1) degenerate Liouville fields and 
a factor u that we absorbed in the definition of 0^-. Finally, we stress that Zq and Zq are 
partition functions of free field theories. They agree with those of the and Liouville 
theory, respectively, if and only if we consider the theory on a surface of genus 3 = 1. For 
higher genus, our relation ( [4.14 ) shows that the partition function Z H of the model is 



related to a 2g — 2 point function in Liouville field theory. 
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4.2 Free field theory computations 

In this subsection, we explain how to compute the residues of the first order poles in WZNW 
correlation function from free field theory (see also our discussion at the end of section 2). 
We then determine the corresponding quantities for correlators in Liouville field theory and 
show that the results agree with our relation (4.14) for the full correlators. 

As we have sketched in section 2.3 it is possible to compute the residues of poles in the 
i?3~ correlation functions by inserting powers of the screening charges into correlators of a 
linear dilaton <j) and a f3"f system. For the rest of this section we shall fix the world-sheet 
metric such that p = 1. The vertex operators and the usual screening charge take the same 
form as above 

Vj{n\z)= l^i+V^V 6 ^ 1 ^ , S= [ d 2 wS{w) = - f d 2 wPPe 2b ^ w '^ . (4.15) 



The TV-point correlation function has a pole at Ylvti" + 1) = 1 — 9 — s (s £ Z>o), whose 
residue is obtained by integrating the following correlators over the positions 

IN s \ N IN s \ 

( n y^Mzu) n s(w k ) ) = n \^ +2 ( n e^+^-^ n \ x 

\u=l fc=l / v=\ \u=l k=l I 

IN s \ I N s \ 

x (n e ^ w n/ 3 ( w fc))(n e "^ w ii[-^*)]) • 

\v=l fe=l / \l/=l k = l I 

(4.16) 

Throughout this entire subsection, correlation functions are properly normalized such that 
the expectation value of the identity is trivial rather than the partition function. Since the 
free boson <p is subject to a background charge Q = = b, the contribution from the 
linear dilaton theory can be computed as 

IN s \ I N s \ 

/ JJ e 2b(j„+l)4»{zu,z„) J-J e 2b^>{w k ,w k ) \ _ ^ goi / j-j- e 2b(j v +l)^ q (zu,Zu) j"J e 2b^> q {w k ,w h ) \ 

\u=l k=l I \v=l k=l I 

(4.17) 

where the twisted zero mode of <p contributes the factor 

N s 
A gol _ TJ gSfcOV+^^solC^,^) TJ e 2b(f> ao i(wk,ui k ) (4-18) 

u=l k=l 

The correlation function of the single valued fluctuation field <j) q = (j)—(j) S oi may be expressed 
through the functions F and H, see appendix |B|, 

IN s \ N s 

/ rj jbUu+i)M*>x>) n e ^H q { Wk ^ k) \ = -q H{Zv) ^{ 3v+ i) -q H ( Wk) ^ x (419) 

\u=l k=l I u=l k=l 

N N s s 

v<\x v =\k=\ k<l 
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The factors in the second line of eq. ( 4.16| ) may be evaluated using the same formulas we 
employed in our path integral derivation. In particular, with the help of formula ( |4.5| ) we 
can conclude that 



N 



\u=l 



k=l 



n^^n^)) = n 

/ k=l 

s 

= n 



N 



k=l 



u=l o=l 



Utilizing the relation between the prime form E and the special function F (see appendix 
IbT) we conclude 



N 



N 



\u=l 



n ^ vi{zv) n ^ ) ( n e-*^ ni-^*)] 

/ \v=l k=l 



2s 



n 

fc=i 



Uu=i\E(z„w k W 



fc=l 



(4.20) 



Equations ( 4.18| ), ( f4.19| ) and ( f4.20| ) provide all the information we need in order to deter- 
mine the residue (|4.16|) of the WZNW correlator at ^2 u (j v + 1) = 1 — 9 — s. 

We would like to rewrite the correlation functions in model in terms of Liouville 
theory. The background charge of Liouville theory is assumed to be Q = 6 + 1/6. We use 
the following vertex operators and Liouville screening charge, 



V a (z,z) = e 2a ^ z ~ z) , S= d 2 wS{w) 



d 2 we 2by{w,w) 



(4.21) 



with a = b(j + 1) + 1/(26). According to the general result (|4.14| ), it should be possible to 
reproduce the residues computed in the previous subsection from the expression 



N N+2(g-l) s \ 

\Q 9 AT \ 2 ( "Pf g2a„v(zf,Zf) "PI" e ~j;<p(.Vi,Vi) e 2btp(w k ,w k ) 



n 

\v=l 



n 



(4.22) 



k=l 



where 0^ is given in eq. ( 4.1C| ) and the correlation function is evaluated in a linear dilaton 
background with background charge Q = = 6+1/6. Using the same formulas as in the 
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previous subsection, we find 



/ N N+2(g-l) s 

YJ ^,a v (p{z v ,z v ) J~J e -\<f>(Vi,Vi) 

\v=l i=l k=l 



N N AT+2( ff -l) 7V+2(g-l) 

i/</i i>=l i=l i<j 

N s N+2(g-l) s s 

w=lfc=l i=l k=l k<l 



N N+2(g-l) s 

x n ^) 2(62+l)(>+l)+l+ ^ n Hivr^u 11 ^ 

v=l i=l k=l 



Our aim now is to replace the factors F{z u ,yiy v+1 which do neither appear in nor in 



the residues of the WZNW model. Using the relation (4.7) with p = 1, we find 



N N+2(g-l) N N 

n n E(z v , yk y» +i tr(*v) 2j '' +2 = n n e ^ ^ +jp+2 > ^ 

u=l k=l v=l fj,<p 

or, equivalently, 

AT N+2(g-l) 

n n n^ yk r +1 H{z v f^ = 

v=\ k=l 

N N 

= {\u\ 2 e s °) s+9 - 1 Yl | Ai |2>+2 e 2fe(>+2)0 so i(^,2.) "Q F ^ Zp y»+J P +2 . (4.24) 
In deriving this equation we have made use of the equality ^2 u (jv + l) = l — g — s. Inserting 



(4.24) into the linear dilaton correlator we rewrite the residues of correlation functions in 
Liouville theory in an appropriate form. Once these are multiplied with our function 
Q 9 N these reproduce the results of the previous subsection and thereby confirm nicely the 
outcome of our general path integral derivation. 

5. Conclusions and Outlook 

In this work we proposed a new and elegant path integral derivation for the correspondence 
between local correlation functions in H% WZNW model and Liouville field theory. Our 
results reproduce the findings of |?J for correlators on the sphere and generalize them 
to surfaces of arbitrary genus g. Correlation functions of the H% WZNW model are 
determined from Liouville field theory through eq. ( 4. 14| ) . Physical correlators are obtained 



in the limit —* of vanishing twist parameters and after integration over w a . A 



correspondence between residues of correlation functions is encoded in eq. ( 4.14 ) and it 



was verified explicitly through free field calculations. For the torus, we have also explained 
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how to map the BPZ equations for Liouville theory to the KZB equations in WZNW 
models. The extension of this analysis to higher genus was not addressed. We believe that 
an explicit comparison along the lines of section 3.2 is possible, though rather cumbersome. 



It seems appropriate to add a few comments on our 2g — 1 moduli w a and to 



explain their relation with the twist parameters introduced by Bernard [23]. We recall 
that Bernard's constructions involve a 3g-dimensional space of twists. It parametrizes a 
set of g group elements which determine boundary conditions of currents along the f3- 
cycles. Our coordinates correspond to very special twists with group elements of the 
form gk ~ exp(— 2iri\k J°) where J° is the Cartan generator. Consequently, insertions of 
the current J°(w) into correlators of the WZNW model may be converted into the action 
of some differential operator. The latter has exactly the same form as in Bernard's work 
(see eqs. (4.16) and (4.17) of [|23}|). For insertions of the component J~(w), the story is a 
bit different. In this case, we may use the relation J~(w) = (3{w) along with our equation 
fl3~3|) to derive 

IN \ H / N 9-1 \ 

\ u=l I (A,w,t) \v=1 <t=1 / 

where flf x ^ r) = / JJ V 3 Mz u )\ . (5.2) 

\ v=l I (A,ro,r) 

On the right hand side, the complex numbers w a that multiply the holomorphic one-forms 
appear in place of differentiation with respect to twist parameters in Bernard's work. In 
this sense, our w a are dual twist parameters. We are not convinced that insertions of 
the third current J + (w) can similarly be replaced by the action of some operator. Let us 
point out, however, that this is not crucial for a successful match between the BPZ and 
KZB-type differential equations. In fact, only certain combinations of the KZB equations 
appear in this context. The relevant ones emerge from inserting the Sugawara tensor at 
the points r/i at which J~ vanishes. Hence, the term J + J~ drops out. 

As explained e.g. in p3] (see also references therein), one of the ramifications of the 
geometric Langlands program involves conformal blocks of WZNW models at the so-called 
critical level and their relation with certain classical W algebras. In the case of the 
model, the critical level is k = 2. Hence, we reach this point in the limit b — > oo in which 
the associated Liouville theory becomes classical. For genus g = correlation functions in 
the if^-Liouville correspondence, the limit of infinite parameter b was analyzed in detail in 
Q. It might be rewarding to carry out a similar investigation for surfaces of higher genus 
g 7^ 0. In this case, the Gaudin Hamiltonians that emerge from the critical WZNW model 
on the sphere get replaced by Hamiltonians of Hitchin's integrable system. 

The relation ( 4.14 ) between correlation functions in the H% model and Liouville field 



theory may be regarded as an 'off-critical' (and non-chiral) version of the geometric Lang- 
lands program. Let us recall that conformal blocks of the WZNW model diagonalize the 
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action of the current algebra on the fusion product of its representation spaces in the same 
sense in which Clebsch-Gordan maps (block-) diagonalize the action of a Lie algebra on 
tensor products. When we are dealing with Lie algebras, the geometric Langlands program 
achieves more: it provides a distinguished basis in the tensor product consisting of eigen- 
vectors of a classical W-algebra. Once we go off-critical, the classical W-algebra becomes 
quantum. In the case considered here, the W-algebra is the Virasoro algebra. As usual, 
the action of the W-algebra is block diagonalized by its conformal blocks. Putting all this 
together, an off-critical version of the geometric Langlands program should single out a 
distinguished basis for WZNW conformal blocks which may be expressed directly through 
conformal blocks of the W-algebra. Our main result eq. ( [4.14 ) claims that for the H% 
model such a basis is given by the WZNW correlators on the left hand side. Let us stress 
that the proper basis is found for (twisted) correlators on any closed Riemann surface. It 
seems likely that a similar off-critical version of the geometric Langlands correspondence 
exists for other Lie algebras (see also comments below). 

There are several extensions of our results that seem worthwhile being analyzed. To 
begin with, it would be interesting to study correlation functions on surfaces with bound- 
aries. For the WZNW model, maximally symmetric boundary conditions were found in 
1 26]. Using new boundary theories for /?7 systems (see [27]) along with ideas from a forth- 
coming paper on branes in the GL(1|1) model [28], a first order formulation for correlators 
with insertions of both bulk and boundary operators can be developed. An evaluation along 
the lines we presented above should then relate these to correlation functions in boundary 
Liouville theory ||, g ||. For some disc amplitudes, such relations between correlators 
on a surface with boundary have been proposed in [32]| . 

More importantly, it is very tempting to address generalizations to WZNW models 
of rank r > 1. First order formulations for models with higher rank are certainly known 
(see e.g. [33] and references therein) and it is likely that these may be employed to reduce 
correlators of WZNW primaries to correlation functions in conformal Toda theories. The 
evaluation of the corresponding WZNW path integral, however, requires significant new 
ideas, mainly because the nilpotent part of higher rank algebras is no longer abelian. This 
is directly linked to a non-linear dependence of the Kac-Wakimoto like action functionals 
on some of the fields 7. We plan to return to these issues in the near future. 
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A. Differential equations for genus one case 



In this appendix we would like to demonstrate that the differential operators ( 3.22 ) and 
( 3.30 ) are related through equation ( 0.33[) . As we described in the main text, it is important 
to first replace the derivatives d^J = d/dz u which are evaluated for fixed \i v in terms of 
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dy. This is achieved with the help of formula ( 3.34 ). In order to verify the latter, we must 
show that 

dm,, = , SiX = . (A.l) 



The second equation can be shown trivially with eq. fl3.9p . The first property of <5j may be 
established as follows. With eq. ( |3.9| ) we find 

lh, X kliv = ^(Vi,Zu)(^2^(zu,yj) - ^£{zv,Zn)) +^2^(yi,Zn)C(zu,z lx )- (A.2) 

3 i^v 
H j+i 

The right hand side is a double periodic function of yi, which could become singular at 
Ui ~ yj, z v . We can analyze the singular behavior of /i" 1 ^/^, with the help of the following 
expansions, 

1 



dzi{z, z) 
t{z,z>? 



z — Z' 



7 -2(z-z') m + 0((z-z') 2 ) , 



1 



(z - z') 2 

1 

{z - z'f " 



-2 Vl + 0(z - z') , 
4r?i + 0(z - z') 



(A.3) 



at z ~ z'. In fact, using the above expansions for £ one can show that \i~^bi[i v has 
no singularities when two of the insertion points yi and approach each other. Since 
the whole expression is double periodic and free of singularities, the function should be 
constant, independent of yi, z v . Therefore, it suffices to calculate it at one single point. Let 
us set yi = z v , then we find 



[z v - Zu) 



'{Zy - yj) 



^ O(zu-z^) j-f. 9{z v - yj ) 



(A.4) 



Here, we have used the following expansion around t ~ z v , 



r(t - e\z v - z„) 



e{t - zu) e(z u - zu) 



+ 



e {z v - Zjj) _ 

9(z u - z^ V 0(z v - Zu) 



{t-z v ) + 0{(t-z v f) 



Furthermore, we can set each of yj (j ^ i) to be one of z^ (/i ^ v) because the above 
quantity does not depend on yj either. From the above equation we can see that the 
equation \i~ x biii v = is indeed satisfied. 

There appears another derivative in the differential operator D|* for torus correlation 
functions, namely the derivative d T with respect to the modulus. We would like to check 
that it acquires no corrections when we switch from the variables fi u to yi,u, or, in other 
words, 

d 



dr 



. 



(A.5) 
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In order to see this, we compute 



^ ^ = E Jl^y + ^ - E -^-^ • (A.6) 
Here we should notice that 

Amd -jw =e ~M = u[z ' o) + 0)2 ' Am& ~iw = ~ 6m + o{z) • (a - t) 

These formulas show that \i~^d T [i v has no singularities. As before we may evaluate the 
derivative at e.g. yi = zi and then find indeed /j, u ~ 1 d T [i u = 0. 

In order to show (|3.33j ), we have to understand the properties of the function @' N that 
we defined in eq. ( 3.32Q , 



N N N 

N 



e' N = e'{o)w H 0(z u - z ^ Y[0( yi - Vj )w H 6{z v - yi yw . (Ai 

u<fi i<j v,i=l 

Conjugation of the various derivatives by the factor O^y gives 

@, *~ l ik & » = ^ + 2^E^^)-2^E^^)> 

d d 1 1 

N d T 9'(0) 1 d T 6{z v - z„) 1 ^ d T 9( yi - yj ) 1 ^ d T e{z u - Vi ) 
2b 2 9>(0) 26 2 £- 0(z„-z„) 26 2 ^- 26 2 ^ 0{z v - Vi ) ' 

ty < ~~.f-l Mil 

Moreover, we note that 



+ 5 E^(^'%)-E^(^'^) E^'%)-E^(^'^) ' ( A - 9 ) 
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and 



Q' N 1 6i®' N 



d 1 



This long list of equations puts us into the position to finally prove eq. ( 3.33 ), 



(evr^QV?) - A L = - ^(£e(w, Vi) -E^.^)) 2 + (a.io) 
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There could be double or single poles at j/, = z u ,yj, but we can check that such terms 
are absent. Moreover, even if we regard the expression ( |A.1Q ) as a function of yj (j ^ i), 
there are no singular terms. Therefore, the problem is whether the constant independent 



of yi, yj vanishes or not. We set yj = Zj for j ^ i so that ( |A.10|) becomes 
{® l NV T 1 Df{Q' N r 1 )-D\ = 

\y% ~ Zi)\ 



1 3 1/ 

{((Vi, Zif + di((yi, Zi)) + ^-f rji - [ 6771 + 
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9(Vi ~ Zi) 



(A.ll) 



where the last term of the right hand side comes from T and in the derivation we have also 
used 

d 

2-iTirj(T)~ 1 —rj(T) = -r]i . 
Taking the limit of yi — » Zj, we indeed obtain the desired equation ( |3.33| ). 
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B. Theta functions on a general Riemann surface 



In this appendix we summarize some basic results concerning a free boson on a Riemann 
surface of genus g. Basically, our exposition follows the discussions in [24]. See also 

\EMM- 



B.l The prime form 

We consider correlation functions on a compact Riemann surface E of genus g with a 
complex structure. Let us choose a canonical basis of homology cycles a k , f3 k (k = 1, • • ■ , g) 
satisfying 

ui = hi , <b ui= m , (B.l) 

where = 1, ■ ■ ■ ,g) denote the holomorphic one- forms on S. The complex symmetric 
matrix r k [ (Imr^ > 0) is known as the period matrix. We now fix an arbitrary point po in 
£ and construct a map from the universal cover E of the surface £ to C 9 , 

z k {p) = [ P u k (B.2) 

with p,po on E. This embedding of E into C 9 is known as the Abel map. Functions on the 
image of the Abel map descend to the surface E if they are periodic under all shifts of the 
form z' k = z k + m k + Tkin l with integer coefficients m k , n l . In order to construct a few basic 
objects on the surface E and its cover E, we recall the following definition of Riemann's 
theta function, 

O s (z\t) = ^ViA^ + Sifmin + dx) 1 + 2{n + 5 l ) k {z + 5 2 )k] , (B.3) 

neZB 

where 8f. = (5ifc,#2fc) with Sik,S2k = 0, 1/2 denotes the so-called spin structure along the 
afc and cycles. Under shifts along the 2g fundamental cycles, 6$ behaves as 

6s(z + rn + m|r) = ex.p[-i7r(n k T k in l + 2n k z k )] exp[27ri(<5f m k - 5%nk)]0s(z\T) . (B.4) 

The Riemann vanishing theorem asserts that 8(z, r) vanishes in a point z on E if and only 
if there exists g — 1 points pi on E such that z can be written in the form 

z = A-Y,Pk , (B.5) 

k=l 

where A is a fixed divisor on the surface E that is known as Riemann class. The right 
hand side of this equation could be considered as an element of C 9 through application of 
the Abel map. Let us now introduce the following holomorphic 1/2-differential hg 

(hs(z)) 2 = d k 9 s (0\r)uj k (z) . (B.6) 



- 30 - 



hs is the essential building block for the important prime form E 

E ^ - ST^TT (B ' 7) 
hs{z)hs{w) 

which is defined for any odd spin structure 5. The prime form E has weight (—1/2,0) x 
(—1/2,0) and near its unique zero at z = w one finds E(z,w) ~ z — w. Moreover, E is 
periodic under shifts z\ along the a^-cycle as z[ = Zi + rii with n[ = 5^ k . On the other 
hand, a non-trivial phase appears if we shift Zi along the /?fc-cycle as z[ = zi + Tik, 

E(z + r k ,w) = - exp (-inT k k - 2ir J E(z,w) . (B.8) 

On the left hand side of this equation, the objects r k denotes the k th column t. k of the 
period matrix Ti k . 

B.2 Free linear dilaton theory 

Let us employ the prime form and some close relatives thereof to spell out the TV-point 
functions of a free bosonic field with background charge Q. For the fluctuation around the 
zero mode, the correlation functions can be given as 

/n e2<W (^)) = il Fi^zj)- 2 ^ \\H{ Zi ) 2 ^ , (B.9) 

\ i=l I i<j i 



where we assume that 



We have defined F and H as 



= Q(l-g) . (B.10) 



F(z,w) = exp I -2-TrIm I u k (Imr)" 1 Im I J J \E(z, w)\ 2 , (B.ll) 

V J w J w / 

and 

H(z) = \p(z)\exp(^- J d 2 w^n(w)\n(F(z,w))^ . (B.12) 

Integration over the g holomorphic forms uj k furnishes an element of C 9 that can be mul- 
tiplied with (Imr) -1 . As in [Q] we can rewrite the function H in a form 

H(z) = exp ( -^-Im [ A u k (Im-r^Im 1^ J) \o(z)\ 2 (B.13) 
\9-l J{g-i)z ' J(g-l)z J 

which involves the Riemann class A that was introduced in the previous subsection. Let 
now p k denote g arbitrary points on £. Then the function cr(z) satisfies 

a{z) _ 9 (z - Y,Pk + A) y-r E(w,p k ) 

11 E( y. m.\ ' V ■ ) 



a(w) 9 (w - J^Pk + A) A-L E(z,p k 
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It is important to mention that cr(z) is a (7/2-differential and that it has no zeros and poles. 
When translated with the k th column vector of the period matrix, a satisfies 



a(z + T fc , w ) = exp iri(g - l)T kk - 2ni / u k 
\ J(g-l)z J 



a(z,w) . (B.15) 
As before, we combined the matrix elements T\ k for I = 1, . . . , g into an element r k 6 C 9 . 
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